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THEORY OF OPEN-ENDED CIRCULAR CYLINDRICAL MICROWAVE CAVITY 


by Norman C. Wenger 
Lewis Research Center 

SUMMARY < 2,0 ° 

The TEqu mode of oscillation in an open-ended circular cylindrical microwave 
cavity is discussed. The cavity consists of a circular waveguide terminated at each end 
with a thin cylindrical partition coaxial with the circular waveguide. The resonant fre- 
quency of the cavity is computed by using Laplace transform and Wiener-Hopf techniques. 
Numerical values for the resonant frequency are presented. 


INTRODUCTION 

Microwave cavities have been used for a number of years to measure the relative 
dielectric constants of gases (refs. 1 to 3). The relative dielectric constant is given 
by the square of the ratio of the resonant frequency when the cavity is evacuated to the 
resonant frequency when filled with gas. 

In many applications, such as in atmospheric research (ref. 4), it is desirable to 
make dynamic measurements of the dielectric constant by having the gas pass through 
the cavity. This requires replacing the solid end walls of the cavity by a termination 
that will totally reflect the microwave energy and yet present very little obstruction to 
the flow of gas. 

Adey (ref. 5) developed a cylindrical cavity with a high Q in which the termination 
took the form of a ring concentric with the body of the cavity and attached to the body of 
the cavity by three spokes. In this design, 67 percent of the end area was removed. 
Thompson, et al. (ref. 6) found that, by judiciously selecting the diameter and thickness 
of this ter minatin g ring, over 90 percent of the end area can be removed without signifi- 
cantly reducing the cavity Q. 

Thom and Straiton (ref. 7) have shown that an open-ended cavity can be constructed 
by simply terminating a section of waveguide at each end with partitions so that each sub- 
division formed by the partitions is a waveguide operating at a frequency below cutoff. 

The partitioned section acts like a perfect reflector to the microwave energy if the parti- 
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Figure 1. - Open-ended microwave cavity. 

tions are sufficiently long. The cross- 
sectional obstruction of the terminations 
for this type of cavity is a result entirely 
of the finite wall thickness of the partitions 
plus any supporting structure. A cross- 
sectional obstruction of less than 2 percent 
of the end area has been reported (ref. 7). 

The open-ended microwave cavity 
shown in figure 1 was recently developed 
at the Lewis Research Center to measure the density of liquid hydrogen under flow condi- 
tions. The density of liquid hydrogen can be computed from its relative dielectric con- 
stant by using an equation such as the Clausius -Mossotti equation. The cavity consists 
of a circular waveguide terminated at each end with a coaxial cylindrical partition. The 
effect of the cylindrical partition is to separate a portion of the waveguide into a coaxial 
waveguide plus a smaller circular waveguide that also serves as the inner conductor for 
the coaxial waveguide. 

This report presents an analysis of open-ended circular cylindrical cavities of the 
type shown in figure 1 to determine the relation between the resonant frequency and the 
cavity dimensions for the case where the cavity is evacuated. The analysis is restricted 
to the TEq^i resonant mode since this mode of oscillation has a very high Q and there- 
fore is widely used. The solution for the resonant frequency of the TEq.^ mode can be 
easily computed if the reflection coefficient of the TEq^ circular waveguide mode inci- 
dent on the cylindrical partition is known. 

The model that will be used to compute the reflection coefficient is shown in figure 2. 
The various regions of interest have been numbered for ease in reference. The model 
consists of a perfectly conducting circular waveguide of radius b and of infinite extent 
in the z-direction. Coaxial with the waveguide is an infinitely thin, perfectly conducting 
circular waveguide of radius a that extends over the range 0 < z < °°. The reflection 
coefficient will be computed for the case of the TE Q1 mode incident from the left. 



Figure 2. - Cylindrical partition in waveguide. 
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Before proceeding with the analysis, it will be instructive to consider the various 
types of waves that can exist in the different regions. It will only be necessary, however, 
to consider the circularly symmetric TE Qn modes since both the TEqj incident mode 
and the model possess circular symmetry. Thus, the field components of interest will 
be the 0 -component of the electric field and the r- and z-components of the magnetic 
field. All other field components will vanish. 


b 

c 

E 0 (r,z) 

Ej(r,z) 

E|(r,z) 

<?(r, 0) 

* + (r, 0) 
*'(r,0) 
G(r,r 0 ,z- z o) 
£(r, a, 0) 

Sf + (a, a, 0) 

SHa, a, 0) 
H r (r, z) 
H z (r,z) 

J in 

J 0 (a, z) 
j 

/(a, 0) 


SYMBOLS 

radius of circular partition (fig. 2) 
radius of circular waveguide (fig. 2) 
velocity of light 
0 -component of electric field 

0 -component of incident electric field 

0 -component of scattered electric field 
bilateral Laplace transformation of E® 

c 

single -sided Laplace transformation of E q over the interval z > 0 

single-sided Laplace transformation of E® over the interval z < 0 
Green's function 

bilateral Laplace transformation of G(r,a,z) 

Wiener-Hopf factor of gr(a,a, 0) 

Wiener-Hopf factor of Sf(a,a, 0) 
r -component of magnetic field 
z -component of magnetic field 
imaginary part of 

0 -component of electric surface current density on surface r = a 
for z > 0 

bilateral Laplace transformation of Jg 
free -space wave number 
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L length between partitions (fig. 6) 

p(/3) function of /3 that gives & + and gT proper behavior as |/3| - » 
R reflection coefficient of TEq^ mode 
|R| modulus of reflection coefficient 
r radial coordinate (fig. 2) 
r Q dummy coordinate 

die real part of 

t time 

z axial coordinate (fig. 2) 

z Q dummy coordinate 

<* n propagation constant of TE Qn mode in region II 

/3 Laplace transformation variable 

0 n propagation constant of TE Qn mode in region I 

0J real part of 

/3^’ negative imaginary part of /3^ 

r n eigenvalue of TE Qn mode in region I 

y n eigenvalue of TE Qn mode in region n 

6 n eigenvalue of TE Qn mode in region HI 

6 polar angle (fig. 2) 



M 0 magnetic permeability of free space 
P n propagation constant of TE Qn mode in region m 
cp phase of reflection coefficient 

a) angular frequency 


WAVEGUIDE MODES 
Region I 

The general solution for the electromagnetic field in region I is of the form 
(ref. 8, pp. 69 to 72) 
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-j/3-.z jfri z V 0 n z 

E 0 (r,z) = A^I^rJe 1 + RA^iyje 1 + > y^^Je 


(la) 


Z i 

n=2 


H r (r,z) = 


A 

CO/i 


oo 

.0. -i^z RA 1 ^ 1 j/3 1 z \ AS 6 2 

Illj^rJe 1 +^-lj 1 (r 1 r)e 1 +) -SJLj^r r)e n 

n=2 


(lb) 


A i r i 

H (r, z) = - L JnO^rte 
Z o 


-j6^z RA 1 r 1 


jwfi. 



6 z 

J 0 ( r n r ) e n ( lc ) 


where the A fl are complex amplitude constants. A time dependence of e^ is implicit 
in these equations. The propagation constants 6 n must satisfy the equations 


g2 _ ,2 _ p2 

Pi - k o 1 1 


0n = r n “ k n n > 1 

^n n o 


(2a) 

(2b) 


where k Q = u>/c, since the general solution must satisfy Maxwell's equations. The 
eigenvalues T n are determined by the boundary condition that requires Eg(r, z) to 
vanish at the perfectly conducting surface r = b (see fig. 2). Thus, 

J l( r n b > = 0 (3) 

which has the roots Tjb = 3. 8317, I^b = 7. 0156, . . . , r^b « (n + 1/ 4 )ir (ref. 9, 

p. 166). 

The first two terms in the solutions (la) to (lc) are recognized as the incident and 
reflected TE Q1 circular waveguide modes, respectively, where R is the reflection 
coefficient. The remaining terms correspond to the TE Qn circular waveguide modes. 

It has been assumed that k Q is in the range < k Q < so that the TEqj mode is a 
propagating mode and the TE Qn (n > 1) modes are evanescent or cutoff modes. 
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Region II 


The general solution for the electromagnetic field in region n is of the form 
(ref. 8, pp. 69 to 72) 


E 0 (r,z) = 



B n J i(r n r ) e 


H r (r, z) 



H z (r, z) 



(4a) 


(4b) 


(4c) 


where the B n are complex amplitude constants. Since the general solutions must 
satisfy Maxwell’s equations, the a» n are given by 


a 


2 

n 




( 5 ) 


The eigenvalues y n are determined by the boundary condition that requires E 0 (r, z) to 
vanish at r= a (see fig. 2). Thus, 


Jl(r n a) = 0 (6) 

which has the roots y 1 a = 3. 8317, y 2 a = 0156, . . . , y R a « (n + l/4)7r (ref. 9, p. 166). 

The terms in the solution are recognized as the TE Qn circular waveguide modes. 
It has been assumed that k Q is less than y^ so that all the TE Qn modes in region n 
are cutoff modes. 


Region III 

The solution for the electromagnetic field in region III is of the form (ref. 8, pp. 77 

to 80) 
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(7a) 


E 0 (r, z) = 


- y c „[ n i< 6 „ 




a)J 1 (6 n r) - 


H. 


(r,z)= 7—^- [ N i( 6 n a)J l (6 n r) “ J l (6 n a)N l (6 n r) ] e ^ 


L^ im o 

n-1 


(7b) 


H z (r,z)= -V' ^2. [N 1 (6 n a)J 0 (6, 1 r) - Jj( V> N 0<V>] e ^ 


(7c) 


n=l 


where the C n are complex amplitude constants. The condition that the solution must 
satisfy Maxwell’s equations requires 


n 




( 8 ) 


The boundary conditions require E^(r, z) to vanish at r = a and r = b (see fig. 2). The 
boundary condition at r = a is satisfied automatically by equation (7a). The boundary 
condition at r = b requires the eigenvalues 6 n to satisfy the equation 

N l (6 n a)J l (6 n b) ' J l (6 n a)N l (6 n b) = 0 ( 9 ) 

The roots of this equation are dependent on the ratio b/a. For example, if b/a = 1.2, 

6ja= 15.7277, 6g a = 31.4259, ...» 6 n a » mr/Rh/a) - l](ref. 9, p. 204). 

Hie terms in solution (7) correspond to the TE Qn coaxial waveguide modes. It 

has been assumed that k is less than 6. so that all the TE„ modes in region m are 

o 1 on 0 

cutoff modes. 


SOLUTION FOR REFLECTION COEFFICIENT 
Integral Equation for Electric Field 


Since the electric and magnetic fields are related by Maxwell's equations, only the 
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electric field needs to be determined to specify the total field uniquely. Thus, the 
9 -component of the electric field, which is the only nonvanishing component of the 
electric field, can be regarded as a potential function from which the magnetic field 
components can be derived. 

The formal solution for the reflection coefficient R will be obtained by setting up 
an integral equation for the scattered electric field due to the TEq^ mode incident on 
the cylindrical partition. 

Green's function. - The integral equation will be formulated by employing the 
Green's function, G(r, r Q , z-z q ), defined to be the solution of the differential equation 


3G + £g + /j 

Sr 2 r8r dz 2 '° 



6(r - r 

_6(z - z ) 

r 


( 10 ) 


that satisfies the boundary condition 


G(b,r o ,z-z Q ) = 0 (11) 

The Green's function can be thought of as being the 6 -component of the electric field 
radiated from a filamentary loop antenna of radius r Q located at z = z q in an infinitely 
long circular waveguide of radius b. 

The solution for the Green's function will be obtained by using a mode -expansion 
technique. A general solution of equation (10) that satisfies boundary condition (11) is of 
the form 


G(r, r o,z- z o )= D n (r o’ z - z o )J l (r n r > 

n=l 


( 12 ) 


where the D n are expansion coefficients to be determined, and the T n are given by 
equation (3). 

Substituting the general solution into equation (10) gives 



SJ l( r n r ) l 3J l( r n r > <2 l\ 

— + ~ - L ^ L - + K ~ -1 Ji (T n r) 


Sr 


Sr 


+ Jl( r n r) 



n=l 


-6(z - z Q ) 


6 (r - r J 


( 13 ) 
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It can easily be shown that the quantity in brackets is equal to (k 2 - rj^J^r^r). Thus, 
equation (13) can be reduced to the form 


oo 

Z 

n=l 


3 2 D 



(k 2 - I> 

2 

\ 0 iy n 

’bTT 



6(r - r ) 

J l( r n r ) = - 6(2 - 2 q ) — 


(14) 


A differential equation for each of the D n can be obtained by using the orthogonal 
properties of the Bessel functions. Multiplying both sides of equation (14) by rJ^r^r) 
and integrating with respect to r over the interval 0 ^ r ^ b give 


3 2 D. 


n 


Sz 



- T 2 )d 

n / n 


2 J l^ r n r o^ 

1)2 [Jo< r n b) ? 


6(z - z Q ) 


(15) 


The solution to equation (15) must be continuous at z = z Q and have a discontinuity in 
the derivative with respect to z at z = z Q of an amount -(2/b 2 )|jj(r n r o )^[jQ(r n b)] 2 |. 
In addition, the solution must correspond to an outward propagating or decaying wave at 
z = ± 00 . The particular solution that meets all these requirements is 


D (r , z-z ) = 
n v o’ o' 




J l< r n r o> 


\1/ 2 0 

) t>[W>] 


2 




(16) 


/ 2 2\^/ 2 / 2 2\^ 2 

where the branch of (T* - k^j must be selected so that 9te (T^ - k I > 0 and 
/ o 2\^/ ^ 

./mlT - k^j >0 to satisfy the conditions at infinity. If the expression for D r is 
substituted into the general solution (12), the Green's function can be expressed in the 
form 


GK r , r o,z-z o ) 


J l( r i r o) J i( r ir) -j|3 1 (z-z o ) 
e 

tflb 2 [J 0 < r l b >] 2 



J l (r n r o )J l (r n r) e -' 5 n (z ' z o> 

^Vofrnh )] 2 


z > z n 
o 


(17a) 
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J l< r n r o> J l< r n r > A (z ' Z o> 
e z < z Q 

V 2 t J o < r „ b » 2 

(17b) 

As |z - z q | becomes large, the value of G is given by the first term in either 
equation (17a) or (17b). 

Electric field . - The 6 -component of the electric field E^(r,z) must satisfy the 
differential equation 


G( r , r Q , z-z Q ) = J l (r i r o )J l (r i r) e ^l (z ~ Zo) + 

j^bV^b )] 2 



0 + i — d - + — e - + (k 


3r 


2 r 3r 


3z 


2 --L) 

° 2/ 


E 0 = O 


(18) 


and the following boundary conditions (see fig. 2): 


E 0 (b, z) = 0 all z 
E 0 (a, z) = 0 for z > 0 


(19a) 

(19b) 


In addition, the electric field must correspond to the incident TEq^ mode as z — - c 
Thus, 


-j/3i z 

lim E 0 (r,z) = J^l^rje (20) 

The reflected TEq^ mode does not appear in equation (20) since it will later be assumed 
that the medium within the waveguide has small electrical losses so that the reflected 
field will have a negligible amplitude at z = 

Integral equation . - Multiplying equation (10) by E q and equation (18) by G and 
then subtracting give 


e 0 v 2 g 


G = - 6(z - z Q ) 


«(r - r ) 


E 0 (r,z) 


( 21 ) 
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where 



r^o Next, each term in equation (21) will be inte- 

Figure 3 . - Surface s. grated over the volume V bounded by the 

surface S, as shown in figure 3. The various 
portions of the surface have been denoted to Sg for ease in reference. The sur- 
faces Sp Sg, and Sg are at infinity. Performing the integration over the volume V 
gives 


f y (e 0 V 2 G - GV 2 E„)dV= -2.E e (r 0> z 0 ) 


( 22 ) 


The integral over the volume V can be reduced to an integral over the closed bounding 
surface S by using Green’s theorem (ref. 10, pp. 803 and 804) with the result 



-2nE 9 (T 0 ,z Q ) 


(23) 


where n is the outward normal coordinate to the surface S. 

Since both E 0 and G are zero along the perfectly conducting surface r = b 
(boundary conditions eqs. (11) and (19a)) the portion of the surface integral over Sg 
will vanish. The surface integral over Sg and Sg will also vanish since E 0 is zero 
at z = +°° due to all the TE Qn modes in regions II and HI being cutoff modes. In 
addition, the first term in the surface integral must vanish along the surfaces Sg and 
S^ since Eg is zero along the perfectly conducting surface r = a for z > 0 (boundary 
condition eq. (19b)). Thus, the only contributions to the surface integral will come from 
the integration over Sj and from the second term in the integrand over Sg and S^. 

It can easily be shown by using the asymptotic forms of Eg and G as z — -°° 

-j/3.Z 

that the integral over the surface Sj is equal to -2i7Jj(r^r Q )e •. Along the sur- 
faces Sg and S 4 , the quantity dEg/dn is equal to -3E 0 /3r and 3E 0 /3r, respectively. 
Equation (23) can therefore be expressed in the following form: 
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E 0< r o’ z o } = J l( r i r o) e 


-J01 


y»°o 

dEn 

e (r.z) 
3r 

3Eo 

- — (r, z) 
9r 


I G(a,r o ,z-z o ) 


/o 


r=a+ 

r=a- 


dz 


(24) 


In order to simplify the form of equation (24), the function Jg (a, z) will be defined 


as 


J q (a, z) — 



- 


“ 

i 

3Eg(r, z) 

3Eg(r, z) 


jcuM 0 

3r 

3r 




r=a + 

r=a“ 


(25) 


It can be easily shown by using Maxwell's equations that J g(a, z) is the electric surface 
current density on the surface r = a for z > 0. The value of Jg(a, z) is zero, of 
course, for z < 0. Thus, equation (24) can be rewritten in the form 


Eg(r, z) = Jj^jrJe 




- ]wM o a 


r°° 

J G(r,a,z-z o )J 0 (a,z o )dz o 


(26) 


where r has been interchanged with r , and z with z Q for clarity. The symmetry 
property of the Green's function G(a, r, z Q -z) = G(r, a, z-z Q ) (see eqs. (17)) was used to 
obtain equation (26)). Equation (26) shows that the total electric field Eg is the sum of 
the TE Q1 mode incident field given by the first term on the right plus an integral that 
corresponds to the field radiated by the electric current on the surface r = a for z > 0. 

It will be convenient in the following analysis to split the electric field Eg into 
two parts: an incident field Eg and a scattered field Eg so that 

Eg(r, z) = Eg(r, z) + E®(r, z) (27) 

where 

i -j/M 

Eg(r, z) = Jj(r jr)e 1 (28) 
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Combining equations (26) to (28) gives the desired integral equation for the scattered 
electric field E®: 


E®(r, z) = -ja>M o a 



a,z-z o )J 0 (a,z o )dz o 


( 29 ) 


In going from equation (26) to equation (29), the lower limit on the integral was changed 
from 0 to -«o since J 0 (a, z Q ) is zero in the range -°° < z Q < 0. 


Solution of Integral Equation 

Laplace transforms. - The solution of the integral equation for the scattered 
electric field will be obtained by using Laplace transform and Wiener-Hopf techniques. 
Let the functions g(r, a, /3), and / (a, 0) be the bilateral Laplace transforms with 

respect to z of E®(r, z), G(r,a, z), and J 0 (a, z), respectively: 


g (r, 0) = f E«(r, z)e”^ z dz 

OO 

(30a) 

9 (r, a, 0) = f G(r, a, z)e _/3z dz 

•4.00 

(30b) 

/(a,/3)= f Jg(a, z)e"^ z dz 

*4.oo 

(30c) 


In order to make all the Laplace transforms exist in a common region in the com- 
plex /3-plane, the propagation constant will be made complex. This is equivalent 
to introducing an energy loss mechanism into the medium interior to the waveguide. Let 
/3f = /3j - j/Sj’ where > 0 and 0^ >0. In the final solution, /3j will be set equal to 
zero to recover the solution for the lossless case. 

The functions 6 and J are, at this point, unknown since E® and J 0 are un- 
known. However, it can be shown by using the expressions for the general forms of the 
solutions given by equations (1), (4), and (7) that & is analytic in the strip -/3" j 

and that J is analytic in the region 3Pe/ 3 > -0y 

The transform of the Green’s function can be computed at once by taking the trans- 
form of the solution for G given by equations (17). This approach will yield the trans- 
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formed Green's function in the form of an infinite series. It will be more convenient in 
the following analyses, however, to compute the Green's function by first transforming 
the basic equation and boundary condition for the Green's function given by equations (10) 
and (11) and then solve the transformed equations. This approach will yield the solution 
for the transformed Green's function in closed form. 

The Laplace transforms of equations (10) and (11) are 


a 2 s? , l S9 1 

3r 2 r3r 



5(r - a) 
a 


(31a) 


9(b, a, 0) = 0 (31b) 

where the variable r Q has been set equal to a to conform with equation (30b). It should 
be noted that the transform of G is taken with respect to z with z Q = 0. The solution 
of equation (31a) that is bounded in the interval 0 < r < b and satisfies equation (31b) is 
of the form 


9(r, a, 0) = A(a, 0)Jj(Ar) r < a 


(32a) 


Sf(r, a, 0) = B(a, ^fj^AbjN^Ar) - J^ArjN^Ab)] r > a 


(32b) 


2 X / 2 
where A = (k + p ) 
o 

The functions A(a, 0) and B(a, 0) are determined from the conditions imposed by 
the right side of equation (31a), which requires to be continuous at r = a and to have 
a discontinuity in the first derivative with respect to r at r = a of an amount -l/a. 
Solving for A and B and then substituting the results into equations (32) give 


<?(r, a, 0) = - - 
2 


J 1 (Ab)N 1 (Aa) - J 1 (Aa)N 1 (Ab) 


Jj(Ab) 


J i (Ar) 


r ^ a 


(33a) 


9(r, a, 0) = - - 
2 


J 1 (Ab)N 1 (Ar) - J 1 (Ar)N 1 (Ab) 
J i (Ab) 


J i (Aa) 


r > a 


(33b) 


Since 9 is an even function of A, the choice of either branch of A will yield the same 
result for 9, It can be shown that 9 is analytic in the strip -0" <9te0 < 0". 

The Laplace transform of the integral equation for the scattered electric field 
given by equation (29) is 
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£{r, 0) = -ju>M o a^(r, a, 0)/(a, 0) 


(34) 


where the faltung theorem (ref. 10, pp. 464 and 465) has been used to take the transform 
of the integral. It will be convenient to express the transform of the scattered electric 
field g as the sum of two single -sided transforms / + (r, 0) and g~{r, 0) where 


/ + (r, 0) 


*"(r, 0) 




E|(r,z)e"^ z 


dz 


dz 


(35a) 


(35b) 


The function S~ is, at this point, unknown. The function g + , however, can be 
computed at the radius r = a by using the fact that the scattered electric field at r = a 
for z > 0 must simply be the negative of the incident electric field since the total electric 
field must be zero to satisfy the proper boundary condition. Using the expression for the 
incident field given by equation (28) gives 


/ 


£+(a, 0) = - / E^(a, z)e dz = - 


/ 


JjtTjaJe 


£ + i0i 


(36) 


The transform of E^(a, z) is only defined in the region 3fe0 > -0^’. 

Combining the result for /+ with equations (34) and (35) yields a form for the 
transformed integral equation that can be solved by using the Wiener -Hopf technique: 


*-(a,0) 


J^a) 

P + iPi 


ja>M o a3(a, a, 0)/(a,0) 


(37) 


It should be noted that, when the previous equations were combined to obtain equation (37), 
the radius r had to be set equal to a in each of the equations since could only be 
computed at r = a. 

Wiener -Hopf technique . - Equation (37) is a transformed integral equation of the 
Wiener -Hopf type (ref. 10, pp. 978 to 980). This equation can be solved for g~ if the 
transformed Green's function can be decomposed into the factors 9f + and 9~ so that 
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<3 = <S Jr l<g- , where <3 + is analytic and nonzero for Step > -0V and is analytic and 

f f ^ 

nonzero for Step < 0j. This step allows equation (37) to be rewritten in the form 

J (r a)sr( a ,a,/3) 

$ ‘(a, /3)sT(a, a, p) - — ± — t = -jam aS^(a, a, 0)/(a, p) (38) 

13 + 3 /^ 0 

The left side of equation (38) is analytic in the region Step < p'^ except for the pole at 
P = introduced by the second term, whereas the right side of equation (38) is 
analytic for Step > -p'^. The left side can be made analytic everywhere in the region 
Step < p 1 ^ by adding the term [3 {&, a, -iP^)]j(P + 30!) to each side of the equation. 
Thus, 


<r(a,0&-(a,a,0)+ J^a) 


ST (a, a, - j/3 j) -sr(a,a,0) 
0 +30j 


Jj^ja^^a^, -j0 x ) 

0 + 30i 


jo)M o a^ + (a,a, p)/{ a, 0) 


(39) 


Now, the left side is analytic for Step < 0*1 and the right side is analytic for Step > -0i'. 
The equality in equation (39) holds only in the strip -p'^ < Step < 0j . 

A function &\p) can be defined from equation (39) so that &\p) is equal to the left side 
of equation (39) for Step ^ -0’,’, the right side for Step > 0V, and to either side in the strip, 

tt !? 1 * 

-01 <Slep < 0j . Thus, $\p) is analytic everywhere in the finite complex 0-plane. 

It can be shown by using the asymptotic expressions for and for 

large 0, derived in the following section, that ^"(0) must also approach zero as 0 — 
Consequently, #~(p) must be zero everywhere since zero is the only function that is 
analytic everywhere and vanishes at infinity (Liouville’s theorem). Since 3tP) is zero, 
each side of equation (39) must also be zero. Thus, 


<T(a,0) 


Jl(ri a )[sr(a, a, -j0i) - ST (a, a, 0)] 
S?"(a,a,0)(0 + j0 x ) 


(40a) 


/(a, 0) 


J x (r i a )S?“(a, a, -50j) 
jwM o aST f (a, a, 0)(0 + j0j) 


(40b) 


16 



Edge conditions . - The conditions that have been specified for 9* and 9~ are that 
9= 9+/9~, 9 + be analytic and nonzero for #e0 > -0p and 9~ be analytic and nonzero 
for #e0 < jSj*. These conditions do not uniquely specify this pair of functions since both 
9+ and 9~ can be multiplied by any function that is analytic and nonzero everywhere in 
the finite 0 -plane to yield a new pair of functions that will satisfy all the original condi- 
tions. Any such pair of functions when combined with equations (40) will yield a solution 
for the scattered electric field that satisfies all the conditions imposed on the field up to 
this point. 

The remaining boundary condition that must be imposed to specify the field uniquely 
is commonly called the edge condition (ref. 11). The edge condition requires the total 

electric field E fl when evaluated at r = a to be of the order z*/ 2 as z — 0. The edge 

y - 1/2 
condition also requires the current density Jg to be of the order z ' as z — 0. 

Equation (28) shows that the incident electric field Eg is of order 1 as z — 0. 

Thus, the scattered electric field must also be of order 1 as z — 0 since the scattered 
field must have a component of this order to nullify the effect of the incident field so that 
the total field can be of the order 7 }^ as z - 0. The edge conditions can therefore 
be stated as 


E® (a, z) = 0(1) z - 0 (41a) 

Jg(a, z) = 0(z -1 / 2 ) z - 0 (41b) 

It can be easily shown that, if a function is of the order z a as z — 0 (a >-l), its 
Laplace transform with respect to z will be of the order as 0 — °°. Con- 

sequently, the edge conditions require the transformed quantities and f to satisfy 
the conditions 


< ?-(a,0) = <7(0- 1 ) 0-°o 

(42a) 

/(a,0) = tf(0' 1/2 ) 0-co 

(42b) 


Combining equations (40) and (42) reveals that the functions 9 + and 9~ must be of 
the order 0 -1 / 2 and 0 1 / 2 , respectively, as |0| — °°. 

Wiener -Hopf factors of Green* s function . - The transformed Green’s function eval- 
uated at r = a is given by (see eq. (33)) 


Sf(a, a, 0) 


Ji(Aa) 

[j^XtyN^Aa) - J 1 (Aa)N 1 (Xb)] 

2 Jj(Xb) 11 11 


(43) 
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where A = ( k Q + $) • The decomposition of # into the ratio of g? + to is best 

carried out by expressing 9 in the form of an infinite product. 

The infinite product expansion of J^(Xa)/Xa is given by (ref. 10, pp. 382 to 385): 


J^(Aa) 

Aa 




( 44 ) 


since the function J^(Aa)/Aa is equal to 1 at A = 0, has a zero derivative with respect 
to A at A = 0, and has simple zeros at A = ±y n (see eq. (6)). The infinite product 
expansion of J^(Ab)/Ab can be obtained from equation (44) simply by replacing a with 
b and y n with r . Thus, 


Jj_(Ab) 

Ab 



(45) 


The infinite product expansion of the function J^(Ab)N^(Aa) - J^(Aa)Nj(Ab) is given by 


J 1 (Ab)N 1 (Aa) - J 1 (Aa)N 1 (Ab) = - (- 

7 r \a 




(46) 


since this function is equal to 2/ir(b/a. - a/b) at A = 0, has a zero derivative with respect 
to A at A = 0, and has simple zeros at A = ±6 n (see eq. (9)). 

Combining equations (44) to (46) reveals that the transformed Green's function can be 
written in the infinite product form 



Equation (47) can be put into a more useful form by reintroducing the parameters 0 n , 
a n , and p n from equations (2), (5), and (8), respectively. Thus, 
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Equation (48) shows that the transformed Green’s function has simple zeros at 0 = ±a R 
and /3 = ±p n and simple poles at /3 = ±j/3 1 and /3 = ±0 n (n > 1). Each of the infinite 
products can now be easily factored into the product of two infinite products so that one 
infinite product is analytic and nonzero in the region 9le \ 3 < and the other is analytic 

and nonzero in the region . The result of this operation allows and S? - 

to be easily identified as 



(49a) 



2 2 

The quantity 1 - a /b has been arbitrarily associated with the factor S?+. The 
exponential terms in equations (49) were introduced to make the infinite products con- 
verge. The logarithmic derivative of each of the infinite products generates an infinite 
series. If the exponential terms were absent, the n^ 1 term in each series would be of 
the order of 1/n; consequently, each series and the associated infinite products would 
diverge. The exponential terms require the n** 1 term of the series to be of the order of 
1/n , which ensures that both the infinite series and infinite products will converge. The 
form of these exponential terms is not unique since they are only required to generate 


19 



terms of the order 1/n that will exactly cancel the existing terms of order 1/n in the 
infinite series. 

The function p(0) in equations (49) must be selected to give the functions 27+ and 
ST the proper algebraic behavior for large 0, as discussed in the previous section. In 
order to solve for p(0) it will be necessary to determine the asymptotic forms of the 
infinite products for large 0. 

Consider the infinite product 



a n + f Y-fta/ror 
, r n ) 


Since a n and y n approach mr/a as n 


°°, it is apparent that the infinite product 




-0a/mr 


only differs from the previous infinite product by a bounded function of 0. Consequently, 
the asymptotic forms of the previous two infinite products will be identical for large 0. 

The gamma function T(x) can be expressed in infinite product form as (ref. 10, 
pp. 421 and 422) 


OO 



where y is Euler’s constant (y = 0. 5772157. .. ). The asymptotic form of T(x) is 
(ref. 10, p. 424) 

T(x) =* (27T) 1 / 2 x x -( 1 /2) e - x (51) 

as x - °°. Combining equations (50) and (51) shows that 



(i +5) e - x/n ~ (2 tt r 1/2 x -t x+(1/2) le (1 ' r)x 


(52) 
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as x — °o. 


Thus, the initial infinite product 



-jSa/mr 


is of the order 


^y[W*W/2)] e (l-y)/3a/ir 

as /3 — This same procedure can be applied to each of the infinite products appearing 
in equations (49) with the result 


sr^a.a,/}) = exp 


Sktof-k-V-lnf— ' \|p(0 

ir \b - a / it \b - a/ 


as 0 — °°, and 


ST( a,a,P) = 0P 1/ ' 2 




(53a) 


(53b) 


as /3 Since 9 + and 9~ must be of the order of 0 _1//2 and /3 1 / 2 , respectively, as 

| /3 1 — «> it is apparent from equations (53) that a suitable choice for p(/3) is 


p(/3) = exp 



(54) 


Thus, the solution for the Wiener-Hopf factors of the transformed Green’s function is 
completed. 

Scattered electric field. - The scattered electric field evaluated at r = a can now 
be computed by inverting the transformed field £ : 



— / *(a, /3)e^ z d/3 

2aj Jq 


(55) 
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The inversion contour C must be located in the 
strip -j3j 3 < ft^, as shown in figure 4, 
since this is the only common region in the 
/3-plane in which all the transforms are analytic. 

The transformed field £(a, /3) is the sum of 
<? + (a, ft) given by equation (36) and /"(a, ft) given 
by equation (40a). Combining equations (36) and 
(40a) with equation (55) gives 

a, -jftj) ^ 
(/3 + j^)sr (a, a, ft) 

(56) 



To evaluate the field for z < 0, the contour C can be closed in the right half /3-plane 
with a semicircle of infinite radius, as shown in figure 4. It can be easily shown that the 
integral over this semicircle is zero. Thus, the integral over the original contour C 
must equal -2nj times the sum of the residues of the poles of the integrand in the right 
half /3-plane. The integrand has poles in the right half /3 -plane due to the zeros of 
at ft = j/3^ and at ft = ft n (n > 1). Performing the integration gives 


J 1 ( r 1 a )^"(a, a, -j0j) j/3 x z 


2j/3j MI (a, a, ft) 


1 a/3 

0=301 



J 1 (r 1 a)^-(a, a, -j/3j) 


A z 


(0 + j/3,) — (a, a, ft) 
n 1 a/3 


z < 0 


(57) 


n=2 




n 


Comparing equations (la) and (57) shows that the first term in equation (57) corresponds 
to the reflected TEqj mode, and the terms in the summation correspond to the evanes- 
cent TE Qn (n > 1) modes. 

To evaluate the field for z > 0, the contour C can be closed in the left half 
/3-plane with a semicircle of infinite radius, as shown in figure 4. It can be shown that 
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the integral over this semicircle is also zero. Thus, the integral over the original 
contour must equal 2irj times the sum of the residues of the poles of the integrand in 
the left half 0 -plane. 

The integrand in equation (56) has only one pole in this half plane located at 0 = -j0^. 
Thus, 

_ -j0,z 

E® (a, z) = -J 1 (T 1 a)e 1 z > 0 (58) 

Equation (58) shows that the scattered electric field evaluated at r = a for z > 0 is 
simply the negative of the incident field. This result should be expected. Since the inci- 
dent field does not satisfy the proper boundary condition at r = a for z > 0, the 
scattered field must contain a term of the same form as the incident field for z > 0 to 
nullify this improper solution. Terms that correspond to the various TE Qn modes that 
exist in the region z > 0 do not appear in equation (58) since the electric field associated 
with these modes must vanish when evaluated along the perfectly conducting surface 
r = a for z > 0. 

From this point on, the propagation constant will be regarded as being real so 
that the subsequent results will apply for the case where there are no losses. 


Reflection Coefficient 


The total electric field in the region z < 0 can now be obtained by adding the inci- 
dent field given by equation (28) to the scattered field given by equation (57) with the 
result 


-j0iZ 

E 0 (r,z) = J^TjrJe 


ST( a, a, -j0j) 


2101 — (a, a, 0) 
30 


J0lZ 

JjO^rJe 


0=j0! 


3T(a,a,-j0 1 )J 1 (r i a) J^r) 0 n z 


(0 n + J0i) — ( a > a > 0) 

n 1 90 


J l< r n a > 


z < 0 


(59) 


n=2 


0»0„ 


I 
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In obtaining equation (59), the amplitude factor of each TE Qn mode in the scattered 
field was multiplied by J ^ ^ r ' n r )/ ^ 1 ( r n a ) reintroduce the radial dependence of the 

modes. 

If equation (59) is compared with the general solution for the field in the region 
z < 0 given by equation (la), the reflection coefficient R of the TE Q1 mode can be 
easily identified as 


R = 


2?-(a, a, -j/3j) 


2j%^ (a, a, /3) 


1 0/3 

/N01 


( 60 ) 


Combining equation (60) with the expression for <g~ given by equation (49b) allows R to 
be put into the form 


(P n + i/M -j/Bjb/rw 


R= -e 


-2j/3i b A U V n 


a n ' i^l\ 3^i a / n7r 


r, 


n 


/p n ' 3/3 l\ 3/3^(b-a)/njr 


p(i/3i) 


n=l 


n=l 


j^b/rnr 


n + -i^ 3 1 a / n7r 


<P n +3/3 1 \ -j/3 (b-a)/nj7 

n Me 1 PO^) 


6 


n=2 


n 


n=l 


n=l 


(61) 


A careful examination of equation (61) shows that, apart from the first term 
-exp(-2j/3jb/jr), which has unit modulus, the expression for R is simply the ratio of two 
complex numbers with the numerator being the complex conjugate of the denominator. 
Thus, it must follow that |r| = 1. 

This result should be expected. The TEqj mode incident on the coaxial cylindrical 
partition shown in figure 2 (p. 2) must be perfectly reflected since it has been assumed 
that all the TE Qn modes in regions II and III are cutoff modes and that there are no 
losses. Thus, the reflected TEqj mode must have the same amplitude as the incident 
mode. 

The argument or phase of R that will be denoted by cp must equal the sum of twice 
the phase of the numerator in equation (61), since the numerator is the complex con- 
jugate of the denominator, plus jr, which results from the initial minus sign. Thus, 
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Phase of reflection coefficient, <?, rad 



+ 



7T 




( 62 ) 




Figure 5. - Phase of reflection coefficient against dimensionless 
frequency for constant values of radius ratio. 


The terms in the summations are 
recognized as the phases of the 
various infinite products that appear 
in equation (61). The last term 
corresponds to twice the phase of 
p(-j^j). It should be noted that the 
exponential terms introduced to 
make the infinite products converge 
did not appear in the final solutions 
for the magnitude and phase of the 
reflection coefficient. 

The phase of the reflection 
coefficient is not unique since any 
integral multiple of 2ir can be 
added to cp without changing the 
numerical value of the reflection 
coefficient. The values of cp given 
by equation (62) will be in the range 
0 <cp < 27 r. 

Numerical values for cp were 
obtained by using an IBM 7094 
computer. The first 100 terms in 
each of the infinite series were 
retained. It was estimated by using 
an integral technique that truncating 
the series at 100 terms introduces 
an error in cp of less than 0. 01 
radian. The numerical results are 
shown in figure 5 in the form of 
curves of cp against k Q b for typical 
values of the parameter b/ a. The 
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range of k Q b is bounded. The lower limit 
is determined by the condition that k Q b 
must exceed Tjb, which has the constant 
value 3. 8317 so that the TEqj mode in 
region I can propagate. The upper limit is 
determined by the condition that k Q b must 
be less than the smaller of y^b or Sjb. 
This condition ensures that neither the 
TE Q1 mode in region II nor the TEq^ 
mode in region HI will propagate. The values of y^b and 6^b are dependent on the 
parameter b/a. It can be shown that, if b/a is less than 1. 831, y^b < 7.0153 < 6^b, 
and if b/a exceeds 1. 831, 6jb < 7.0153 < y^b. For the case where b/a= 1. 831, the 
TE 01 modes in regions n and HI have identical low-frequency cutoff points, which 
correspond to the values y^b = 6 jb = 7. 0153. 

The numerical results for cp show that at the lower limit of k Q b the value of <p is 
equal to ir for all values of the parameter b/a. Thus, at the lower limit of k Q b the 
cylindrical partition reflects the incident TEq^ mode as if a perfectly conducting sur- 
face spanned the waveguide cross section at z = 0 (see fig. 2, p. 2) since a perfectly 
conducting surface gives a reflection coefficient of unit modulus and argument it. As 
k Q b increases, the incident TEqj mode is still perfectly reflected since the modulus of 
R remains at 1. However, since the phase of the reflection coefficient decreases with 
increasing k Q b, it appears as though the position of the perfectly conducting surface 
moves. Thus, the cylindrical partition acts like a M virtual” perfectly conducting surface 
that spans the entire waveguide cross section whose position is a function of k Q b and the 
parameter b/a. 



z * -L z * 0 

Figure 6. - Model of open-ended cavity. 


RESONANT FREQUENCY OF CAVITY 

An open-ended cavity can be constructed by placing two cylindrical partitions of finite 
length within a circular waveguide as shown in figure 6. The partitions need only be long 
compared with 1/ a ^ and 1/p^ for the formula for the reflection coefficient of the in- 
finitely long partition to apply. To determine the conditions for resonance in the cavity, 
consider the electric field in the region between the partitions to be the sum of TE Q1 
modes propagating in the positive and negative z -directions. Thus, 

-j0i z 30i 2 

Eg(r, z) = AJj(rjr)e 1 + BJ^r^e (63) 
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Cavity length to radius ratio, L/b 


Radius ratio, 
b/a 


Ll\U L5 




Radius ratio, 
b/a 

— L831- 


4.5 5.0 5.5 6.0 6.5 

Dimensionless resonant frequency, k Q b 


7.0 7.5 


Figure 7. - Ratio of cavity length to radius against dimensionless reso- 
nant frequency for constant values of radius ratio. 



~ 10 [ 


:Data point for cavity, 
shown in figure 1 

H b— h— I 

Radius ratio, 
b/a 

^ 2 . 082-1 


4.0 4.5 5.0 5.5 6.0 6.! 

Dimensionless resonant frequency, k Q b 

Figure 8. - Ratio of cavity length to radius 
against dimensionless resonant frequency 
for radius ratio of 2. 082. 


At z = O, the ratio of the reflected wave BJ^r^r) to the incident wave AJj(r^r) must 

equal the reflection coefficient R. Likewise, at z = -L, the ratio of the reflected wave 
j0,L -J/3.L 

AJj(rjr)e to the incident wave BJ^r^rJe must also equal R. Since R is 

equal to e^ these conditions are simply 


B=e^ 


-j^L 


= e i<P 


It can be easily shown that, if there is to be a nontrivial solution for the amplitudes A 
and B, <p and L must be related by the equation 





(65) 


j^L 


e 2 H"^ L 


= 0 


which has the solution 


/3jL = <p (66) 

The required spacing between the partitions to resonate the TEq^ mode can be 
found by solving equation (66) for L by using the values of cp given in figure 5. The 
spacing for the TEq-^ resonant mode can be obtained simply by adding (n - l)2ir to cp 
before solving for L. 

The numerical results for the TEq^ mode are shown in figures 7 and 8 in the 
form of curves of L/b against k Q b for typical values of the parameter b/a. The 
range of frequencies over which the cavity can be timed by varying L/b is a strong 
function of b/a. The value b/a = 1.831 gives the maximum tuning range. 

The numerical results shown in figure 8 are for b/a = 2.082. This value is com- 
monly used in practice since it corresponds to locating the cylindrical partition at the 
radius where the electric field of the TEqj mode has its maximum intensity (i. e. , 
J^(Tjr) has a maximum at r = b/2. 082). This value of b/a was used in the design of 
the cavity shown in figure 1 (p. 2) where the radius a was taken as the mean radius of 
the cylindrical partition. A data point corresponding to the measured values of L/b and 
k Q b for this cavity, as shown in figure 8, indicates that the theoretical and experimental 
results are in good agreement. 


CONCLUSIONS 

A theory has been developed for a class of open-ended circular cylindrical micro- 
wave cavities by using Laplace transform and Wiener-Hopf techniques. Numerical re- 
sults were presented for the resonant frequency of the TE q11 mode. The theory was 
shown to be in good agreement with experimental data. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April, 6, 1966. 
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Page 18, 
Page 18: 

Page 18, 
Page 18: 

Page 18, 

Page 19, 
Page 19: 


equations (44), (45), and (47): The right side should be multiplied by 1/2. 

Line 5 should read "since the function Jj(xa)/xa is equal to 1/2 at x = 0, 
has a zero derivative with respect." 


equation (46): The coefficient of the right side should be - 1/ 7r instead of 2/n. 


Line 12 should read "since this function is equal to - 1/ n(b / a - a/b) at X = 0, 
has a zero derivative with respect. " 


equation (47): 


The quantity 



should be 



equations (48) and (49a): The right side should be multiplied by -1/2. 

2 2 

Line 9 should read "The quantity -1/2(1 - a /b ) has been arbitrarily associ- 
ated with the factor t S + . The. " 


Page 20: 


Line 9 should read "Since a n and y n approach (n + 1/4) n/a. as n it is 
apparent that the infinite product." 


Page 20, line 10: 


The quantity 



should be 



(n + 1/4 )tt_ 


Page 20: Equation (52) should be followed by 

"as x — oo. It can also be shown that 



^ x -(x+e+l/2) e (l-r)x 


as x— °° where e is an arbitrary constant. " 

Page 21: Line 1 should read "Thus, the initial infinite product." 

Page 21, line 4: The first exponent should be -[(j3a/jr) + (3/4)] instead of -[(0a /n) + (1/23- 
Page 24, equation (61): The last quantity in the numerator should be p(-j/3j) instead of 
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Page 18, 
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Page 19: 

Page 20: 

Page 20, 


equations (44), (45), and (47): The right side should be multiplied by 1/2. 

Line 5 should read "since the function Jj(Xa)/xa is equal to 1/2 at X = 0, 
has a zero derivative with respect. " 


equation (46): The coefficient of the right side should be -1 /tt instead of 2/n. 


Line 12 should read "since this function is equal to -l/jr(b/a - a/b) at X = 0, 
has a zero derivative with respect." 


equation (47): 


The quantity 



should be 



equations (48) and (49a): The right side should be multiplied by -1/2. 

2 2 

Line 9 should read "The quantity -1/2(1 - a /b ) has been arbitrarily associ- 
ated with the factor 5T 1 ". The." 


Line 9 should read "Since a? n and y n approach (n + 1/4) 7r/a as n — °o, it is 
apparent that the infinite product." 


line 10: 


The quantity 



should be 



j3a 

(n + 1/4 )tt . 




Page 21: Line 1 should read "Thus, the initial infinite product." 

Page 21, line 4: The first exponent should be -[(j3a/jr) + (3/4)] instead of - [(0a/ n) + (1/2)]- 
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( A 

Page 18, equation (47): The quantity [1 - ¥■ — | should be [1 - 
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Page 19, equations (48) and (49a): The right side should be multiplied by -1/2. 

n p 

Page 19: Line 9 should read "The quantity -1/2(1 - a /b ) has been arbitrarily associ- 
ated with the factor y + . The." 

Page 20: Line 9 should read "Since a n and y n approach (n + 1/4) 7r/a as n — it is 
apparent that the infinite product." 


fl ± \ should be 

1 + p* 

V m 

(n + l/4)jr _ 


Page 20: Equation (52) should be followed by 

"as x — «5. It can also be shown that 

1 + — - — \ e -x / n = ^ x -(x+e+l/2) e (l-y)x 
n + e/ 

as x-°° where e is an arbitrary constant. ” 

Page 21: Line 1 should read "Thus, the initial infinite product." 

Page 21, line 4: The first exponent should be -[(/3a/ n) + (3/4)] instead of -[(/3a/jr) + (1/2)]- 

Page 24, equation (61): The last quantity in the numerator should be p(-j/3j) instead of 
pUJSj). 
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Page 18, 
Page 18: 

Page 18, 
Page 18: 

Page 18, 

Page 19, 
Page 19: 

Page 20: 

Page 20, 


equations (44), (45), and (47): The right side should be multiplied by 1/2. 

Line 5 should read ” since the function J^(xa)/xa is equal to 1/2 at \ = 0, 
has a zero derivative with respect.” 

equation (46): The coefficient of the right side should be -1/n instead of 2/n. 

Line 12 should read ”since this function is equal to -l/jr(b/a - a/b) at \ = 0, 
has a zero derivative with respect.” 

2 \ / 2 
equation (47): The quantity (l - £ — | should be |l - — — 


equations (48) and (49a): The right side should be multiplied by -1/2. 

2 o 

Line 9 should read ”The quantity -1/2(1 - a /b ) has been arbitrarily associ- 
ated with the factor # + . The." 


Line 9 should read ” Since Q! n and y n approach (n + l/4)7r/a as n — «>, it is 
apparent that the infinite product.” 


line 10: 


The quantity 



should be 



(n + 1/4 )n _ 


Page 20: Equation (52) should be followed by 

n as x — °°. It can also be shown that 



^x'( x+e+1 / 2 V 1-r ^ X 


as x— °o where e is an arbitrary constant. " 

Page 21: Line 1 should read ’’Thus, the initial infinite product. " 

Page 21, line 4: The first exponent should be -[(/3a/ 77) + (3/4)] instead of - [(/3a/ w) + (1/2)]- 

Page 24, equation (61): The last quantity in the numerator should be p(-jj3j) instead of 
P(j^ 1 )- 
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Page 18, 
Page 18: 

Page 18, 
Page 18: 

Page 18, 

Page 19, 
Page 19: 

Page 20: 

Page 20, 


equations (44), (45), and (47): The right side should be multiplied by 1/2. 

Line 5 should read “since the function Jj(xa)/xa is equal to 1/2 at x = 0, 
has a zero derivative with respect.” 

equation (46): The coefficient of the right side should be -1/n instead of 2/ n. 


Line 12 should read “since this function is equal to -l/?r(b/a - a/b) at X = 0, 
has a zero derivative with respect. ” 

. 2 \ / x 2 

equation (47): The quantity [l - £ — | should be [1 - 


equations (48) and (49a): The right side should be multiplied by -1/2. 

9 9 

Line 9 should read “The quantity -1/2(1 - a /b ) has been arbitrarily associ- 
ated with the factor 5? + . The.” 


Line 9 should read “Since a n and y n approach (n + l/4)7r/a as n — °°, it is 
apparent that the infinite product.” 


line 10: 


The quantity 



should be 



0a 

(n + 1/4) tt 



Page 20: Equation (52) should be followed by 

”as x — oo. It can also be shown that 



^ x -(x+e+l/2) e (l-y)x 


as x -*• °° where e is an arbitrary constant. " 

Page 21: Line 1 should read ’’Thus, the initial infinite product.” 

Page 21, line 4: The first exponent should be - C(/3a/ 7r) + (3/4)] instead of - [(/3a/ tt ) + (1/23- 
Page 24, equation (61): The last quantity in the numerator should be p(-j/3j) instead of 
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